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SUMMARY 

Recent  proposals  by  Oomory  and  others  for  solving  linear 
programs  Involving  Integer-valued  variables  appear  sufficiently 
promising  that  It  Is  worthwhile  to  systematically  review  and 
classify  problems  that  can  be  reduced  to  this  class  and  thereby 
solved.  Hlstorlcalljfc  non— linear,  non-convex  and  combinatorial 
problems  are  areas  where  classical  mathematics  almost  always 
falls.  It  Is  therefore  significant  that  the  reduction  can  be 
made  for  problems  Involving  multiple  dichotomies  and  k-fold 
alternatives  which  Include  problems  with  discrete  variables, 
non-llnear  separable  minimizing  functions,  conditional  con¬ 
straints,  global  minimum  of  general  concave  fxanctlons  and 
combinatorial  problems  such  as  the  fixed  charge  problem, 
traveling  salesman  problem,  orthogonal  latln  square  problems, 
and  map  coloring  problems. 
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ON  THE  SIONIFICANCE  OF  SOLVING  LINEAR  PROGRAMMING  PROBLEMS 

WITH  SOME  I^f^EaER  VARIABLES 

George  B.  D.^ntzlg 

Recently  R.  Oomory  developed  a  theory  of  automatically 
generating  "cutting  planes"  which  permits  efficient  solution 
of  linear-  programs  In  Integers  In  a  finite  number  of  steps 
[  1  ]  .  This  approach  has  been  generalized  to  a  cat*  where  eone 
variables  are  continuous  and  some  are  constr-alned  to  be  Integers, 
t  E.M.L.  Beale  [2],  and  In  a  more  direct  way  by  Gomcry  [j],  see 
also  [4].  Small  scale  test  problems  have  been  successfully 
computed.  The  procedure  Is  so  promising  that  It  Is  I’elevant 
to  systematically  review  and  classify  problems  that  can  be  reduced 
to  this  class.  We  shall  show  that  a  liost  of  difficult,  Indeed 
seemingly  Impossible  problems  of  a  non— linear,  non— convex,  and 
combinatorial  character  ere  now  open  for-  direct  attack. 

The  cutting  plane  appr-oach  was  first  proposed  and  Its 
power  demonstrated  by  successfully  solving  an  example  of  a 
large  scale  tr-avellng  salesman  problem  by  PXilkerson,  Johnson, 
and  th«  author  [3]  •  Manna  and  Markowitz  axplorad  this  taohnlqua 
further  In  [8]  end  pointed  how  It  could  be  applied  to  solve 
problems  Involving  non— linear  objective  forms  (separ-able  In 
the  variables  but  not  convex). 

In  Section  I  we  shall  give  a  general  description  of  the 
cutting  plane  approach  and  then  describe  the  principles  for 
solving  several  general  type  problems.  In  the  later  sections 
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these  will  le  applied  to  several  well-known  prollems.  The 
outline  for  the  paper  is  as  fellows: 

5iectlon  I:  aenerel  Irinc'.  pies 

(a)  The  Ketl.od. 

(b)  Dichotomies. 

(c)  k— fold  Alternatives. 

(d)  Grlecticn  from  man}  pairs  of  regions. 

(e)  Discrete  Variable  Problems. 

(f)  Non— Linear  Objective  irollems. 

(tS)  Conditional  Constraints. 

(h)  Finding  a  jloial  minimum  of  a  concave  function. 
Gee t Ion  II:  Fixed  Ch  r,^e  Problem. 

Gectlon  lII:  The  'I  raveling  Sal«81lftfl  Pr«bl*«. 

.lection  IV:  The  Ci  tho^onal— La  t  In  'quar>e  Prctlem. 

Gectlon  V:  Fou r-CoiorTn-}  0  Map  (If  possible). 


I;  General  1  Inclples 

(a)  Ti.e  McKiod:  The  cutt’n^;  plane  niethon  consists  In 
first  solvin'}  the  llneai  pro}:8mmLnj  problem  without  the 
integer  constraints.  If  the  optimum  solution  h’  ppens  to 
satisfy  these  conditions  all  Is  ^ell.  If  not  then  additional 
linear  inequality  constraints  (called  cutting  planai)  are  added  to 
the  system  in  sucii  a  wa^  as  to  removt;  the  non— edmi  ssl  I  le  extreme 
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point  sol.itlon  and  yet  retain  all  admleelble  solutions  (e.g., 
those  having  Integer  values).  In  principle,  that  this  could  be 
dane  haa  been  known  for  aoaetlBe.  For  exaaple,  a  plane  that  goes 
through  all  nelghbox^ng  vertices  of  the  non— admissible  extreme 
point  can  be  used  as  cutting  plane  In  the  case  where  all 
variables  must  have  values  0  or  1.  However  such  a  procedure 
has  been  regarded  as  probably  too  slow  and  actual  problems 
until  recently  were  solved  using  more  efficient  cutting  planes 
whose  validity  depended  on  pedal  arguments.  Tills  weakness 
has  been  overcome  by  the  recent  proposals  which  generate  cutting 
planes  In  an  efficient  manner.  It  Is  the  author's  belief  that  new  It 
Is  only  a  matter  of  time  before  a  subroutine  for  Integer  and  partial 
Integer  solutions  will  be  part  of  electronic  con^niter  simplex 
codes. 

Let  us  now  turn  to  the  main  subject  of  this  paper,  types 
of  problems  that  are  reducible  to  linear  programs  some  or  all 
of  whose  variables  are  Integer-valued. 

Quite  often  papers  will  appear  In  the  literature  which 
formulate  a  problem  In  L.P.  (linear  programming)  form  e> cept 
for  certain  side  conditions  like  ^^*^2  *  ^  terms 

of  this  type  such  as  x^*X2  ^  which  Imply  for  nennegatlve 

variables  that  at  leaat  one  variable  ef  eaeh  pair  must  be  sere. 
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Superficially  this  seems  to  place  the  protlem  In  the  area  of 
quadratic  pro-^rammln.^.  However  the  presence  of  such  conditions 
can  entirely  change  the  character  of  the  problem  (as  we  shall 
see  In  a  moment)  and  should  serve  a  wsmln::  to  those  who  would 
apply  wllly-nllly  a  general  non— linear  programming  method.  If 
we  graph  the  conditions  =  0,  Xj^  >  0,  x^  >  0,  x-^  +  x^  >  1, 

the  double  lines  depict  the  domain  of  feasible  solutions.  It 
will  le  noted  that  It  has  two  disconnected  parts.  If  there  are 

■ 

I  \ 

- - 

I 

man>  such  dichotomies  in  a  larger  problem  It  can  result  in  a 
domain  of  feasible  solutions  with  many  disconnected  parts  cr 
connected  non— convex  regions.  For  example,  k  pairs  of  variables 
In  which  one  Is  *ero  might  lead  to  2  disconnected  parts.  Usual 
mathematical  approaches  can  guarantee  at  Lest  a  local  optimum 
solution  to  such  problem^  l.e.^a  solution  which  Is  optimum  only 
over  some  connected  convex  pert. 

It  has  teen  well  known  that  by  special  devices  that  the 
local  optimum  solutions  could  be  avoided  In  many  cases  by  the 
Introduction  of  Integer  valued  variables  but  this  has  only  been 
of  passing  Interest  until  the  recent  developments  rendered 
this  approach  practical.  Our  purpose  here  will  be  to 


P-1486 

^-8-58 


this  knowledg* 

(b)  Dichotomies:  Let  as  be^in  with  the  important  class  of 
problems  that  have  "ei the]>or"  conditions,  ror  such  a  problem 
to  be  difficult  computationally'  there  must  le  many  sets  of  such 
conditions.  Let  us  focus  our  attention  on  one  of  them,  say 


(1)  EITHER  a(x^,x^....,x^)  >  0 

(?)  OR  H(Xj,X2,.--,X^)  >  0 


must  hold  for  values  of  (x, ,x„, • • • ,x_ )  chosen  over  some  set  S. 

1  n 

We  do  not  exclude  the  case  of  both  holding  if  possible.  For 
example* a  contractor  in  a  bid  mi^ht  stipulate  either  x^  >  $10, OCX) 
or  -  0.  If  all  bids  are  nonnetsatlv#  so  that  x^^  >  0,  then  we 
can  write 


EITHER  x^  -  10,000  >  0 

OR  “  ^1  >  0  . 

From  other  considerations  it  may  be  known  that  no  tld  can  exceed 
$1,000,000  so  that  the  set  S  of  Interest  is  L  <  x^  <  1,000,000. 

We  now  assume  that  lower  bounds  for  the  functions  3  and 
H  are  known  for  all  values  of  ,X2#  •  •  •  ,x^ )  in  3.  If  Is 
a  lower  bound  for  G  and  for  H  then  for  ■  1  the  condition 

(5)  C}(x^,X2,  •  •  •  ,x^)  -  6L^  >  0 

holds  for  all  values  of  Slmllsrly  for 
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a  0  the  condition 

(4)  H(x^,X2,  •  •  •  )  -  (1  -  6)L^  >  0 

holds  for  all  values  of  (x,  ,x^,,  •  •  •  ,x^ )  In  3.  For  our  example 

X  £.  IJ 

we  would  have 

-  10,000  -  c-)(-10,000)  >  0 

_x^  -  (1  -  (-1,000, 000)  >  0  . 

The  either— or  condition  (1,2  )  can  now  be  replaced  uy 


(5) 

Q{x^,x^,  *  “  *^n  ^  "  “^^G  ^  ^ 

(S  -  0,1) 

(C) 

H(x^,X2,  •  •  •  )  -  (1  -  ^)Ljj  >  C 

(7) 

0  <  8  <  1 

where  ^  Is  an  Inte^jer  variable.  The  effect  of  ^  =  1  la  to  relax 
the  G  condition  when  H  holds  and  of  r**  »  0  la  to  relax  H  when  G 
holds.  If  G  and  H  are  linear  functions  we  have  reduced  the 
el  then— or  condition  to  three  alraultaneous  linear  Inequalities 
In  wh/ ch  the  variable  ^  must  le  0  or  1. 

A  dichotomy  can  be  used  to  describe  an  L-shaped  region 
(non-convex ) :  for  example,  x^  >  ^1  -  ^2  ^ 

either  <  1  or  x^  <  1 .  We  replace  this  by 
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0  <  Xj  <  1  +  6 

0  <  x^  <  2  -  6 

/  • 

0  ^  8  <  1  (8  -  Opl) 

/ 

V  / . 
/  ^ 

If  now  a  problem  contalna  not  one  but  several  such  pairs 
of  dichotomies  (l)  and  (2)«each  one  would  be  replaced  by  a 
simultaneous  set  (5),  (6),  (7)  in  integer  variables  6^^. 

(o)  lUifold  Alternailves;  More  generally  suppose  we 
have  a  sdt  of  conditions 

(8)  0^(*l,X2,...,x„)  >  0 

aj(xj,Xj,...,x„)  >  0 

jXgJ  •  •  •  ^  0  » 

Suppose  a  solution  is  required  in  which  at  laaat  Ic  of  tho 

eendltions  louat  hold  almiltanooualy.  Wo  roplaoo  this  by 

(9)  G^(x)  -  4  >  0 

Gjx)  -  Og  ^2  ^  ^ 

V*)  -»p  I-p  >  0 

where  Lj^  le  the  lower  i.ound  for  (l(x)  for  x  -  (Xj^,Xg,-**,x  ) 
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In  S  and  are  inteper-valued  variables  satisfying 


(10) 

(Sj  +  ^2  ^  S 

(11) 

0  ^  1  1- 

An  example  of  this  type  of  problem  might  occur  If 
one  wishes  to  find  the  minimum  over  the  shaded  regions 
described  by  0^  0,  ^2^0,  0^  ^  0  and  at  least  two  of 

the  conditions  ^  0«  ^  0«  0^  ^  0  as  In  (12). 


(1?) 


/ 

/ 

/ 

\  / 


\ 

\ 


o 
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/  >  ° 
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aj  ..  Oj(x) 
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(d)  Selection  from  many  paira  of  re(^iona;  The  six-pointed  "Stir 
of  David"  region  shown  on  the  right  in  (12)  can  Lest  be  described 
by  a  dichotomy  in  which  a  point  must  be  taken  from  one  of  two 
triangles.  It  is  only  when  there  are  many  such  pairs  to  be 
chosen  at  the  same  time  that  the  prollem  t ecomes  significant. 

In  general  we  might  have  several  pairs  of  rOf^lons 
( i  •  •  *  j  ,  Rjl^ )  and  the  solution  point  x  must  lie  In  either 
or  r|  for  each  i.  For  each  pair  H  and  F'  we  proceed  as 
follows.  Let  region  R  be  described  ty  a  set  of  inequalities 
Q^ix)  >  0,  a2(x)  >  0,.-.,a^(x)  ^  0  and  H*  by  H^(x)  >  0, 

H^{x)  >  C,«*»,H^(x)  >  0.  The  condition  that  the  point  must  be 
selected  from  either  the  first  or  second  region  can  be  written 

(13)  0^(x)  -  >  0  Hj(x)  -  (1  -  S)L^  >  0 

Ggix)  -  ^L2  >  0  H2(x)  -  (1  -  01^2  ^  ^ 

H„(x)  -  (1  -  Ol;  >  0 
0<''<1,  (^*0orl) 

Where  ire  lower  bound!  for  and  The  more  general 

case  of  selection  from  several  regions  can  be  done  by  Intro¬ 
ducing  several  6^  as  In  (lO)  and  (ll). 

(a)  Discrete  Variable  Problems:  Suppose  that  a  variable  Is 
constrained  to  take  one  of  several  values:  x^^  ■  a^^  or  x^  -  a^,..,, 

*1  “  *k  several  other  variables  areas  are 

also  constrained  the  same  way.  It  would  be  a  formidable  task  to 
test  all  the  combinations.  Instead  we  replace  each  k-folJ  dichotomy  by 


(14) 
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Xi  »  01^  +  ®2^?  +  •••  + 

(15)  ■»■  4  -  1  bj  =»  0  or  1. 

Similarly  let  x  =■  ‘  *^n  ^  repre^c  )t3  a  vector  which  may 

only  take  on  specified  vector  values  x  °  a^  or  x  »  a""  or 

X  =»  •  •  • .  This  may  Le  replaced  ij 

(16)  X  =  a^b^  4  a^^>^  4  •••  4  a^^j^ 

(17)  ^1  "*■  '^2  *“  ^  it  “  ^  -  C  or  1. 


This  device  permits  the  replacement  of  a  non-llnear  function 
^IJ  “  *  svstem  ^  “  1,2, ...,m) 

ty  a  sprinkling  of  representative  values  of  x^,  say  Xj  -  Xj 
where  r  »  1,2, •••,k.  In  this  case  the  vector  la  the  set  of 

IT' 

values  (^1  j » i some  value  Xj  *  Xj. 


(f)  Non— Linear  Objective  Ptollems:  Suppose  the  objective 
form  can  be  written 


(1^) 


n 

^  tj(Xj  )  =>  z(Mln  ; 


where  Is  non— linear  and  non— convex.  Let  each  ^;(xj  te 
approximated  a  broken  line  function.  I’hese  define  a  set 
of  Intervals  1  »  l,^;,***,k  of  width  h^^  and  slopes  s^  for  the 
approximating  chords.  We  now  define  as  the  amount  of  over 
lap  of  the  Interval  from  0  to  x  with  interval  1.  Ttien 


P-.1486 

9-S-58 

-II- 


(19)  X  -  -»•  y^^ 

and  ^(x)  iB  2lven  approxlniately  \  y 

(?C)  <)(x)  *  4  Sj^yj^  4  B2y2  Sj^yj^ 

where 


(21)  C  <  ^1  i  “  l,2,--..k. 

In  the  case  of  convex  4^,  the  procedure  Is  to  replace  x  and 
0(x)  by  (19)  and  (20)  and  condltlona  (21).  Here  the  slopes 
are  monotonlcally  Increaaln.y  so  that 

(2  )  »!  <  ®2  <  •••  <  *k- 

For  a  fixed  x,  (>(x)  would  be  minimum  If  Is  chosen  maximum; 
then  (.ylven  maximum,  so  that  Hiaxlmum;  etc.  In  other  words 
for  the  mlnlralzln,^  solution  the  are  the  overlap  of  the  1^^‘ 
Interval  with  the  Interval  C  to  x  and  all  Is  well. 

Ho. -ever  If  ^(x)  Is  not  convex  as  In  (27),  then  simple 
replacement  of  x  and  <t)(x)  would  result  for  fixed  x  In  y^^  with 
smaller  slopes  leln^  maximized  first.  In  this  case  the  segments 
that  comprise  y^^  would  te  disconnected  and  our  approximation 
for  ^(x)  would  no  Ioniser  be  valid.  In  order  to  avoid  this  we 
Impose  the  condition  that 
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(.3) 

EITHER 

‘‘i  -  ” 

') 

OR 

^1+1  ■ 

f 

^ , 

which  Implies 

that  unless 

>1  is  maximum 

that 

>'1^1  -  C  and 

Yj  Is  maximum 

then  Yi^^i  > 

0  Is  possible. 

We 

I'ewrlte  this 

condition 

(i-H) 

EITHER 

^1  -  ^1  i 

0 

OR 

-  yi4i  ^ 

0 

and  then  replace  It  formally  by 


{<b) 


y^—  hf  —  (-hi)  c^i>^  l-l,2,»**,k  —  1. 

-H.l  -  -  i)>0 

■ ,  =  0  or  1 


upon  substitution  of  -  1  —  6^  simplifies  to 


^1+1  -  ^1+1  ®1 
0  <  6'  <  1 


(6^  -  0,1) 


Ihe  aiove  procedure  for  the  non— convex  case  was  discussed  In 
the  paper  of  Manne  and  Markowitz  [  6  ].  The  convex  case  will 
ce  found  in  [?]  and  [8], 
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h,-, 


h, 


0(x  j 
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A  second  method  Is  worth  notln.j  Lssed  on  (I6).  Any 
point  on  the  curve  ^(x)  can  be  represented  as  a  wel:::;hted 
average  of  two  successive  breakpoints.  Hence  we  may  replace 
X  and  ^(x)  by 


(  j  X 
()(x) 
1 


X  a  4- 

X,  a .  +  •  •  • 

+  X,  a, 

00 

1  1 

k  k 

X  b  + 

X,  b,  +  ... 

+  X 1  b , 

0  0 

1  1 

k  k 

+  +  X,  +  •  •  •  4-  X, 
O  1  1 


(0  <  x^  <  1) 


and  then  Impose  the  conditions  that  all  x^  ■  0  except  for  one 
pair  X^,  ^1^1*  k  ■  4  this  may  te  expressed  by 


( ) 


X  <6 
o  i 


h<^0*  ®1 


X^  £ 
< 

^4  1 


+  6^  + 


+  6^  6. 


-♦■5^  +  6^ 

♦  6., 
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where  are  Inte  ^ei- valued  variables  ssilsf^ln,. 

(3c)  “^3  ^ 

Indeed  It  will  ue  noted  that  when  h,  =1  for'  some  1  =  1 

I  c 

o 

that  the  Irequalliies  Involvin'  A,  and  are  relaxed  but 

o  ^  c 

the  r'cmalnder  satisfy  X,  <  0  since  their  =  0  by  (3C). 

(  ,)  Conditional  Constraints;  Suppose  x  and  y  are  functions 

of  several  varlailes  (x,  ,x,  )  for  -whlcli  upper  bounds 

i.  <-  n 

U  and  lower  bounds  L  and  L,  are  known.  We  wish  to  Impose 

X  j/ 

cond 1 tlcns  such  as 

(31 )  X  >  C  y  > 

We  can  write  this  as 

(32)  KrTHhb  X  >  C,  y  >  C 

Ch  X  <  C 

which  we  rewrite  as 

(33  ;  ^ 

>  b  s 

X  <  (1  -  (^  =  C,l; 

wher-e  the  first  Inequality  is  wrl  ‘  --0  (>  1' stead  of  (>; 

Lecause  “^^he  condition  y  >  C  Is  a utcme ■  1  ce  1 1  y  relaxed  for 
X  =  0  selectin'  1. 

We  can  now  elaroi'Ste  this  to  Impose  conditions  such  as 
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(?»)  X  >  0  =»  u  ^  0 

X  <  0  ® 

which  may  be  written  ea 

(35)  »  ^ 

*  ^  ^2  °x 

w  2  Sg 

+  &2  ■  ^  (^1  ■ 

For  example,  suppose  In  a  T-perlod  program  we  wish  to 
complete  a  specified  work  load  by  the  earliest  period  possible. 

^  w 

Let  be  the  cuamulatlve  sum  of  activity  levels  from  the  t 
period  thru  the  last  period  T,  then  we  wish  to  arrange  matters 
so  that  x^  •  0  for  the  saalleat  t.  In  this  case  we  can  define 
for  t  -  1 ,  2,  . . . ,  T 

(36)  ®  *t  "  ® 

which  «e  nay  rewrite 

(37)  x^  i  -  0,1 

where  Is  an  upper  bound  for  x^. 
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and  determine  Min  z  ^  er« 

(  5b )  E  - 

Find  in  a  lot  el  minimum  cf  a  concave  Tunc- lor.*.  Suppose 
the  concave  function  Z  ■  Z(xj  . . .  ,x^ )  Is  to  be  mlnlrrlzed  over 

a  re  ^1  on  I  .  We  shall  assu-te  '  convex  for’  convenience  liere 
noMn  Itrat  the  devices  d’sc.ssed  earlier  extend  Llie  domain  to 
iF.e  wide  class  express  1  ;,lo  ry  elLr.er— cr  conditions.  We  suppose 
;■  to  ;e  ilven  af’er  st^ltarle  criSn  e  In  varbp.les  In  slander'd 
linear  programming  1’cit, 

(  59)  -  x  =  e  ^  ^ 

wh.ere  r.  Is  a  iven  n  x  n  matrix  and  e  a  ^Iven  m— component  vector. 

This  Is  Intrl  nilcellt  a  dlfTlcrlt  pr'oilem  recause  the 
concave  function  could  have  local  minima  at  many,  Indeed  at 
ell  the  extreme  points  of  •. 

IT.e  concave  iTincilcn  -  ma^,  le  Iven  explicitly  oi' 
ie  Iven  implicitly.  For  exa  ’.pl e,  s  uppose  vector'  y  and  quantity 
z  for  f 1 xed  X  I  3  I ven  . : 

(40;  Fy  -  f  I'X  >  2l 

z  =  ax  —  y.ln  Ty 
:<  ix 

vji.e’.f'  t  and  b  are  ,1  ven  matrices  and  f,  u,  ei;d  ,iven  vec..cr3. 

Vr  ;  3  is  tie  sit  rat  loti  'IIS':  ssed  in  're  appllca'icn  of  these 
* ;  r  ’  3  appllca-lcn  le'.elcpeo  Jcln'T  wl  tr.  Prlllp  Wolfe. 
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methods  to  Solving  Two  Move  Games  with  Perfect  Information 
[9]  •  Here,  however,  we  shall  suppose  that  Z  can  reasonably 
be  approximated  at  all  points  x  In  R  by  the  minimum  Z  of  a 
finite  set  of  k  tangent  hyper-planes. 

(41)  Z  -  “ii*!  ^  “12*2  "  -  ’’l  ^  . - 

to  the  surface  Z  «  Z(x).  The  problem  reduces  to  choosing  Min  Z 
where  Z  must  satisfy  at  least  one  of  the  conditions 

(42)  Z  -  +  ajgXj  +  •••  +  a^„x^  -  b^]  ^  0 

Z  -  [a2ix^  +  aggXg  ^  +  ag^x^  -  b^]  ^  0 


2  -  *  »k2^2  *  •••  \n*n  “  ^  0 

which  we  may  rewrite  as 

(43)  Z  -  [a^jx^  ^  a^jXj  +  •••  +  a^„xj  2-«  61  (1-1,2.... 

6^  +  62*'’**  +  &j^*1  (6^=0  or 

where  -  M  Is  some  assumed  lower  bound  for  the  differences;  this 
solution  depends  on  the  approximation  by  k  hyper-planes  of 
the  function  Z  -  Z(x).  The  solution  given  In  [9],  for  the 
case  where  Z  Is  given  Implicitly  by  (4o),  requires  finding  x,  y, 
Min  z,  and  auxiliary  variables  v  =  (v^,T2,  . . .  ,Tj^)  and  11  j  ^  0 
for  J  -  l,2,...,n'  satisfying 


r-:  >!  y'. 

►IS- 


(JU) 


Sx  =  e ,  Py  «  f  Ex ,  z  »  ax  —  Py 
*  ’'j  - 


1 ,2 , . . .  ,n 


EITHEP  I  nj  <  0 


OP 


y^  <  0 


^  H 

where  tt  »  » '^2*  *  *  *  *\' ^  ^  vector,  Is  the  J  column  of 


t  h 

P,  Pj  the  J  component  of 


FIXED  CHANGE  PROBLEM 

Earlier  we  described  a  problem  where  a  bidder  required 
that  either  the  order  x  ■  0  or  x  >  a.  In  this  and  many  other 
problems  there  Is  an  underlying  notion  of  a  fixed  oharge  that 
Is  Independent  of  the  size  of  the  order.  In  this  case  x  -  a 
represents  the  break-even  point  to  the  bidder.  In  general  the 
cost  C  Is  characterized  by 

kx  +  b  If  X  >  0 

(1)  ^  - 

[  0  if  X  =  0 

where  b  Is  the  fixed  charge.  We  may  write  this  In  the  form 
(2)  ^  -  kx  +  6t  (6  -  0,1) 

where  x  =  0  If  b  =»  0  which  we  Impose  by 


(3) 


X  <  6U 
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(4)  0  <  S  1  (^  -  0*1) 

whtr€  U  la  boim  upper  bound  for  x.  A  dlacusalon  of  the  fixed 
chaise  problem  Including  this  device  will  be  found  In  the  paper 
by  Warren  Hlrach  and  the  author  [lO]  . 

III.  THE  TRAVELING  SALESMAN  PROBLEM 

We  shall  give  two  fonoulatlons  of  this  well-known  problem. 
Let  -  1  or  0  according  to  whether  the  t  directed  arc  on 

the  route  la  from  node  1  to  node  J  or  not.  The  oondltione 


(1) 

t  ■  1, . . .  ,n 

(2) 

1  <*  1 , . . .  ,n 

(3) 

tfi  ' 

J  “  l,2,.,.,n 

(4) 

z  ^IJ^IJ  ■  *  (Min) 

express  that  (l)  there  Is  only  one  t  directed  arc,  (2)  there 
Is  one  directed  arc  leaving  node  1,  (5)  there  Is  only  one 

directed  arc  Into  node  J,  (4)  the  length  of  the  tour  Is 
minimum.  It  Is  not  difficult  to  see  that  an  Integer  solution 
to  this  system  Is  a  tour. 

In  the  paper  by  Fulkerson,  Johnson  and  the  author  the  case 
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of  a  symmttrlc  distance  -  d^^  was  formulated  with  only  two 
Indices.  Here  ■  1  or  0  according  to  whether  the  route 

from  1  to  j  from  J  to  1  was  traversed  at  some  time  on  a 


route  or  not. 

The  conditions 

(5) 

2  "ij  ■  2  . 

i 

J  *  1,2,. ..,n 

(6) 

2  “  z{Min) 

express  the  condition  that  the  sum  of  the  number  of 
entries  and  departures  from  each  node  Is  two.  Tliese  conditlona 
are  not  enough  to  characterize  a  tour  even  though  the  x^j  are 
restricted  to  be  integers  in  the  Interval, 

(7)  0  <  Xj^j  <  1, 

since  sub— tours  like 


(8) 


also  satisfy  the  conditions.  However  If  so-called  loop  con¬ 
ditions  discussed  In  [3]  lilce 


(9) 


‘12  ^23  ”*■  ^31  -  ^ 


I 
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art  ImpoMd  (in  tht  samt  manntr  that  cutting  planes  are 
introduced  as  required)  these  will  rule  out  integer  solutions 
which  are  not  admissible. 


IV.  THE  ORTHOGONAL  LATIN  SQUARE  PROBLEM 

A  latin  square  consists  of  n  sets  of  n  objects  (l), 
(2)».**#(n)  asilgned  to  a  n  x  n  square  array  so  that  no  object 
is  repeated  in  any  row  or  column.  Two  latin  squares  are 
orthogonal  such  as 


(1) 

(1) 

(2) 

(3) 

(2) 

(3) 

(1) 

(2) 

(3) 

(1) 

(1) 

(2) 

(3) 

(3) 

(1) 

(2) 

(3) 

(1) 

(2) 

if  the  n  pairs  of  corresponding  entries  are  all  different. 

It  was  conjectured  by  Euler  that  there  are  no  orthogonal  latin 
squares  for  certain  n.  In  spite  of  a  great  deal  of  research 
by  top  notch  mathematicians  the  case  for  n  •  10,  for  example, 
has  never  been  settled.  It  has  been  suggested  informally  by 
David  Oale  that  the  proposed  method  be  tried  in  this  area. 

The  fomulation  is  straightforward  and  well  known.  Let 
*ljkl  “  ^  ^  according  to  whether  the  pair  (i,j)  is  assigned 

to  row  k  colusm  I  or  not.  The  condition  that  the  pair  is  assigned 
to  o.Tly  one  location  is  given  by 
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uTi  ^ 


1>J  *  l(2f>.i>n. 


The  condition  that  at  least  one  Is  assigned  to  each  location 


k,  t  Is: 


^  ^IJk/  ■ 


The  conditions  that  1,  J  appear  only  once  In  the  first  and 
second  latln  square  respectively  In  column  i  Is  given  by 


t  ti  “  1^2f>a»fn> 


J»k  =  lf2>>«>>n> 


Similarly,  the  conditions  that  1  and  J  appear  only  once  In  the 
first  and  second  latln  square  respectively  In  row  k  Is  given  by 


“  1,2, ...,n. 


1,2, •  •  • ,nt 


It  is  Interesting  tc  note  that  every  pair  of  subscripts  that 
are  possible  out  of  four  are  summed  to  form  the  six  sets  of 
n  equations  each.  For  n  10  there  are  600  equations,  which 


are  too  many  for  a  general  linear  programming  code  to  handle 
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at  the  present  time.  However  with  some  short  cuts  Introduced 
It  might  be  tractable  In  the  near  future. 


V.  POUR-COLOR INO  A  KAP  (if  possible) 

A  famous  unsolved  problem  Is  to  prove  or  disprove  that 
any  map  In  the  plane  can  be  colored  using  at  most  four  colors 
where  no  two  regions  that  have  a  boundary  In  comon  (except 
a  point)  have  the  same  color.  We  shall  give  two  ways  to 
constructively  color  a  particular  map  If  possible.  This 
does  not  contribute  anything  to  a  proof  of  the  truth  or 
falsity  of  the  conjecture  except  that  an  efficient  way 
for  solving  particular  problems  on  an  electronic  computer  may 
provide  a  counter  example. 

Without  difficulty  It  can  be  arranged  (as  below)  so 


\ 

Yellow 


/ 

Red 


Black  ■ 


that  three  regions  have  at  most  one  point  In  common  which 
will  be  called  a  node.  ITiere  will  be.  accordingly,  three 
directed  arcs  leading  from  any  node  1  to  other  nodes  J.  It 
Is  well  known  that  If  It  possible  to  four-color  a  map  then 
(and  this  will  be  true  conversely)  It  Is  possible  to  treat 
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the  nodts  as  cities  and  the  arcs  as  routes  between  cities  and 
either  be  able  to  make  a  tour  of  all  the  cities  or 
t#  Bake  a  group  of  mutually  exclusive  sub— tours  of  the  cities 
In  several  even  (sub— cycle)  loops  as  below. 

We  may  associate  with  each  such  even  cycle  sub— tour, 
directed  arcs  that  reverse  their  direction  as  we  pass  from 
node  to  node 


^  ^ 

This  means  the  nodes  1  can  be  classified  Into  two  classes: 
those  which  have  two  arcs  pointing  away  from  them  and  those 
that  have  two  arcs  pointing  towards  them.  Let  us  set  -  1 
am  arc  Is  part  of  a  sub— tour  In  the  direction  of  the  arrow; 
otherwise  x^j  »  0.  Hence 

(1)  0  <  <  1. 

It  Is  understood  that  only  arcs  (l,j)  and  variables  are 
considered  corresponding  to  regions  that  have  a  boundary  In 
common.  All  arcs  (l,j)  that  do  not  correspond  to  boundaries 
are  omitted  In  the  constraints. 


The  conditions 
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(2)  Z  ’'ij  -  26^  (6^  -  0,1) 

J 

express  the  fact  there  must  be  two  arcs  on  some  sub— tour 
leading  away  from  node  1  If  6^  -  7, otherwise  there  are  none. 

The  conditions 

(?)  E  XlJ  ■  2  -  2»1 

1 

state  there  must  be  two  arcs  on  some  sub— tour  leading  Into 
node  1  If  6 -  0,  otherwise  none.  The  three  sets  of  conditions 
(1),  (2),  (3)  are  those  of  a  bounded  transportation  problem 
and  will  be  integers  (at  an  extreme  point)  If  arc  Integers, 
This  would  seem  to  Imply  that  It  Is  only  necessary  to  assume 
that  5^  are  Integers  and  the  x^^  will  come  out  automatically 
Integral  In  an  extremlzlng  solution  without  further  assumptions. 
However  since  the  objective  form  Is  open  to  choice  by  choosing 
It  In  a  non-degenerate  way  It  Is  clear  that  the  extreme  point 
solution  with  integral  would  be  determined  by  the  process. 

A  second  formulatlrn  suggested  Informally  by  R.  Oomory  Is 
straightforward.  Let  the  regions  be  r  »  1,2,...,R  and  let 
tj,  be  an  Integer-valued  variable  such  that 

0  <  <  3, 
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the  four  valuta  t  =»  0,1, 2, 3  corrf'spondlng  to  the  four  colors. 
If  regions  r  and  s  have  a  boundary  In  common  their  colors 
must  be  different.  Hence  for  each  such  pair 

(**)  ‘r  -  ‘s  °  • 

This  nay  be  written  In  elther^r  form: 

(5)  EITIiER  t^  -  t.  >  1 

OR  >  1 

s  r  ^ 


which  we  may  rewrite 


t 

r 


‘o  >  1 


46 


rs 


t 


8 


>  -  5  +  45 


rs 


5 

rs 


(6) 


0,1 


\ 


M486 
REV.  4-9-59 
-C7- 


REPERENCES 

1.  Gomory,  RAlph,  ’'Esttntlalt  of  an  Algorithm  for  Inttgtr 

Solutions  to  Llnosr  Progrtas,  *  *  Bull,  A— r.  Hath.  8oo>« 
vol.  64,  no.  5»  1958* 

2.  Baals,  E.  M.  L. ,  ''A  Mathod  of  Solving  Linaar  Progi'anning 

Problams  whan  soma  but  net  all  of  tha  Varlablas  must 
taka  Intagral  Valuas,*'  Statistical  Taohniquas  Rasaarch 
Group ,  Prlncaton,  approximate  data  ^rch  19^8. 

3.  Gomory,  Ralph,  Sama  subjact  as  Rafaranca  2.  Papar  in 

praparation  based  on  research  at  Tha  RAND  Corporation, 

July  1958. 

4.  Dantslg,  Gaorga  B. ,  On  Integer  and  Partial  Integer  Linaar 

Prograimlng  Problams,  ina  hand  Corporation,  Paper  p— i4i0, 
June  ^0,  i95b. 

5.  DantElg,  George  P. ,  D.  R.  Pulkereon,  and  Selmer  Johnson, 

''Solution  of  a  Large-Scale  Traveling-Salesman  ProMem,  '  ' 
Journal  of  the  Operations  Research  Society  of  America, 
Wovamber  1957,’"vol7T;“no. « ,  p.555-410.'  - 

MarkowltE,  H.  M.,  and  A.  S.  Manna,  ''On  tha  Solution  of 
Discrete  Programming  Problems,*'  Econometrics ,  vol.  25, 
no.  1,  Jariutry  1957,  p.  19. 

7.  DantElg,  Gaorga  P.,  ''Recant  Advances  In  Linaar  Programming.' 

F^agement  Science,  vol.  2,  no.  2,  January  195^#  p.  131^144 

8.  Chamas,  Abraham,  and  Carlton  S.  lemka,  ' 'MinimlEatlon  of 

Nonlinear  Separable  Convex  Functionals,''  Naval  Rasaarch 
^|istics  Quarterly,  vol.  1,  no.  4,  December  19^4,  p.  501— 

9.  DantElg,  George  B.  ,  Solving  Two-Move  Ga^s  with  Perfect 

Information,  * '  Tha  rand  Corporation,  Paper  i*-i459» 

August  11,  1958 . 

10.  Hlrsch,  Warren  M.,  and  George  B.  DantElg,  The  Pixed  Cha^a 

Problem,  The  RAND  Corporation,  Papar  P— 648,  December  iT 


